LU-TP 03-42 
hep-ph /yymmnn 

Study of the Linked Dipole Chain Model 



in heavy quark production at the Tevatron 

A.V. Lipatov 

Physical Department, M.V. Lomonosov Moscow State University, 119992 Moscow, Russi 



,0 



L. Lonnblad 

^ ■ Dept. of Theoretical Physics, Solvegatan 14 A, S-223 62 Lund, Swede^ 

O ' 
<N 

_ ' N.P. Zotov 

^ ■ D. V. Skobeltsyn Institute of Nuclear Physics, 

. M.V. Lomonosov Moscow State University, 119992 Moscow, RussiJ^ 

Abstract 

> 

C ■ We present calculations of charm and beauty production at Tevatron within the framework of kj- 
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factorization, using the unintegrated gluon distributions as obtained from the Linked Dipole Chain 
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model. The analysis covers transverse momentum and rapidity distributions and the azimuthal 
correlations between b and b quarks (or rather muons from their decay) which are powerful tests 



Oh, for the different unintegrated gluon distributions. We compare the theoretical results with recent 

Qh' 

' experimental data taken by D0 and CDF collaborations at the Tevatron Run I and II. 
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I. INTRODUCTION 



It is known that in the description of a given cross section in lepton-proton or proton- 
proton interactions at high energies it is not enough to consider only the leading order 
perturbative terms. Although at large scales, /i, the running coupling constant a s may 
be small, each power of a s is accompanied by large logarithms due to the large phase 
space available for the additional gluon radiation. The solution of this problem is to resum 
the leading logarithmic behavior of the cross section to all orders, thus rearranging the 
perturbative expansion into a more rapidly converging series. 

The most familiar resummation strategy is based on Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) 0, 0, [J |J| evolution equation, where large logarithmic terms proportional 
a" ln n (/i 2 /AQ CD ) are taken into account. The cross sections can be rewritten in terms of 
pro cess- dependent hard matrix elements convoluted with universal parton density functions 
which are described by the DGLAP equation. In this way the dominant contributions come 
from diagrams where parton emissions in initial state are strongly ordered in virtuality. 
This is called collinear factorization, as the strong ordering means that the virtuality of the 
parton entering the hard scattering matrix elements can be neglected compared to the large 
scale fi. 

DGLAP evolution describes most experimental results from electron-proton and proton- 
proton colliders. By using input parton densities which are sufficiently singular when x — > 
0, this formalism can also account for the strong rise of F2 at small x, as observed at 
HERA. However, there are problems with the description of non-inclusive observables such 
as forward jet production in ep and heavy quark production in ep and pp collisions. 

At the energies of the HERA, Tevatron and LHC colliders, the hard scale /i of the heavy 
quark and quarkonium production processes is large compare to the Aqcd parameter but 
on the other hand /j, is much less than the total center-of-mass energy: Aqcd ^ ^ <S s/s- 
Therefore in such case it was expected that the DGLAP evolution should break down. The 
situation is classified as "semihard". 

It is believed that at asymptotically large energies (very small x oc fi 2 /s) the theoretically 
correct description is given by the Balitsky-Fadin-Kuraev-Lipatov(BFKL)P, 0, ?| evolution 
equation. Here large terms proportional to a™ln n (l/x) are taken into account. The BFKL 
evolution equation actually predicts a strong power-like rise of F 2 at small x. Just as for 
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DGLAP, it is possible to factorize an observable into a convolution of process-dependent hard 
matrix elements with universal parton distributions. But as the virtualities (and transverse 
momenta) of the propagating gluons are no longer ordered, the matrix elements have to 
be taken off-shell and the convolution made also over transverse momentum with the un- 
integrated gluon distribution k^). The unintegrated gluon distribution determines the 
probability to find a gluon carrying the longitudinal momentum fraction x and transverse 
momentum kx- This generalized factorization is called " A/p-factorization" IslO]- If the terms 
proportional to a™ \n n (fi 2 / Aq CD ) and a" ln n (/i 2 /AQ CD ) ln n (l/x) are also resummed, then the 
unintegrated gluon distribution function depends also on the probing scale /i. 

The unintegrated gluon distribution in the proton is a subject of intensive studies [lo|. 
Knowledge of these is in particular necessary for the description of heavy quark production 
at future colliders, such as the LHC. This quantity depends on more degrees of freedom than 
the usual collinear parton density, and is therefore less constrained by the experimental data. 
Various approaches to model the unintegrated gluon distribution have been proposed. One 
such approach, valid for both small and large x, has been developed by Ciafaloni, Catani, 
Fiorani and Marchesini, and is known as the CCFM model^, Q, Q]. It introduces 
angular ordering of emissions to correctly treat gluon coherence effects. In the limit of asymp- 
totic energies, it is almost equivalent to BFKL 3, 16, Q], but also similar to the DGLAP 
evolution for large x and high /i 2 . The resulting unintegrated gluon distribution functions 
depend on two scales, the additional scale q being a variable related to the maximum angle 
allowed in the emission. 

The Linked Dipole Chain Model (LDC) jisl . Iioj l is a reformulation and generalization of 
the CCFM model, and agrees with CCFM to leading double logarithmic accuracy. Also 
LDC is formulated in terms of /c^-factorization and unintegrated distribution functions, but 
here these distributions are essential single-scale dependent quantities. In this paper we will 
apply the LDC formalism for the description of the charm and bottom production processes 
at Tevatron. 



The application of /c^-factorization supplemented with the B 



20 



21 



TvL anc 



22 



23. 



CCFM evolution 
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equations to heavy quark hadroproduction is discussed in 
was shown that bb production cross section at Tevatron can be consistently described using 
the /^-factorization approach together with different BFKL or CCFM-like unintegrated 
gluon distributions. The NLO pQCD calculations based on the DGLAP evolution scheme 
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underestimate data by a factor about 2 - 5|2£ 
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321 ] . In a previous paper [27[ the 



dependences of the b quark, B meson and their decay muon cross sections on different forms 
of the unintegrated gluon distribution was investigated. It was found that the properties of 
different unintegrated gluon distributions strongly influences the bb or muon-muon azimuthal 
correlations since these quantities are sensitive to the relative contributions of different 



production mechanisms to the total cross section 
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Based on the above mentioned results here we will use the /cy-factorization approach 



together with LDC unintegrated gluon distribution functions 



experimental data|2 
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]. To illuminate the effect of the different contributions 
we will study three different versions of LDC unintegrated gluon distributions presented 
in[33]. It is interesting to note that the LDC unintegrated gluon distributions has been 
fitted to the inclusive F2 data at HERA, and for heavy quark hadroproduction at Tevatron 
we will obtain essentially parameter-free theoretical results. We also present our predictions 
for the differential cross section of D*, D + and D° meson production in pp collisions at 
Tevatron Run Ilj^. 

This article is organized as follows. In Section |H] we give a short review of the CCFM 
and LDC formalism. In Section 
compare them with the D0 
we give some conclusions. 
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In] we present the numerical results of our calculations and 
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341 ] data. Finally, in Section |^ 



II. THE CCFM EVOLUTION AND LDC MODEL 



In figure ^ we show a typical fan diagram for the initial-state radiation in a DIS lepton- 
proton event. According to the CCFM formalism the emission of gluons during the initial 
state cascade is only allowed in an angular-ordered region of phase space. In addition 
the initial state emissions are ordered in the positive (along incoming proton) light-cone 
momentum k + . All other kinematically allowed emissions (symbolized by q[ emission in 
figure HJ) are defined as final-state emissions. 



The CCFM evolution equation can be written in a differential form 

^ 2 d $(x,k^,q 2 ) _ f dz d(t)P g (z,(q/z) 2 ,k 2 
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dq 2 A s {q 2 ,Q 2 



2tt A s (f,Q 2 ] 



^^(x',k 2 T ',(q/z) 



(1) 



where $(x, k\, q 2 ) is the unintegrated gluon distribution, k T = k' T + cfr, and q = qx/i^ — z) 



lepton 




proton 



FIG. 1: A fan diagram for a DIS event. The quasi-real partons from the initial-state radiation are 
denoted qi, and the virtual propagators k{. The dashed lines denote final-state radiation. 



is the rescaled transverse momentum of an emitted gluon (with azimuthal angle <ft). The 
Sudakov form factor As is given by: 



q 2 1-Qo/q 

a(q 2 {l - z) 2 ) 



— z 



(2) 



with a = 3as/ir. For inclusive quantities at leading logarithmic order, the Sudakov form 
factor cancels the 1/(1 — z) collinear singularity of the splitting function. The splitting 
function P g for branching i is given by: 



P ( z . a 2 k 2 )- ~ Zi)2) + ^3A A ( z . 2 ,2 \ 

1 Z{ Z{ 

where the non-Sudakov (or non-eikonal) form factor A ns is defined as: 

i 

f dz' f dq 2 

lnA ns (z i ,g 4 2 ,/c| i ) = -J — J —a9(k Ti -q)6(q-z'qi). 



(3) 



(4) 



The CCFM equation incorporates both BFKL and DGLAP evolution. In the small x limit 
the unintegrated gluon distribution obeying CCFM evolution equation (0) can be written 
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as: 



where the splitting parameter z is defined as z% = k+^/k+^-i. The interval for z variables 
is between and 1, which guarantees ordering in k + , and the angular ordering condition is 
satisfied by the constraint 

Qi > qi-iZi-t, (6) 

explicitly written out in eq. (J3J). The distribution function depends on two separate scales, 
the transverse momentum kx of the interacting gluon, and q, which determines an angle 
beyond which there is no (quasi-) real parton in the chain of initial-state radiation. It may 
that the role of this variable is similar to that of /i 2 in the collinear gluon 



be argued 
density. 



The LDC model 
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19| relies on the observation that the non-Sudakov form factor in 
equation (j3J) can be interpreted as a kind of Sudakov giving the no-emission probability in 
the region of integration. An additional constraint on the initial-state radiation is added 
requiring the transverse momentum of the emitted gluons to be above the smaller of the 
transverse momenta of the connecting propagating gluons: 

q Ti > mm(k T i, k T i-i)- (7) 

Emissions failing this cut will instead be treated as final-state emissions and need to be 
resummed in order not to change the total cross section. The remaining initial-state gluons 
are ordered both in q + and g_, which implies that they are also ordered in angle or rapidity. 

One single chain in the LDC model corresponds to a set of CCFM chains. It turns out 
that when one considers the contributions from all chains of this set, with their corresponding 
non-eikonal form factors, they add up to one[l||- Thus, the non-eikonal form factors do not 
appear explicitly in LDC, resulting in a simpler form for the unintegrated gluon distribution 
function: 



x 9(q +i ^i - q +i )9(q-i - - JJzi)5(lnfcf. - Ink 



(8) 



Tn) 
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The notation in (8) and what will follow refers to that of figure El Here, a typical DIS 
event is shown together with the corresponding phase space available in the 7*p rest frame, 
where the rapidity y and the transverse momentum qj- of any final-state parton are limited 
by a triangular region in the (y, lngf^-plane. The proton direction is towards the right end 
of the triangle, and the photon direction is towards the left. The real emitted gluons are 
represented by points in this diagram. The virtual propagators do not have well defined 
rapidities, and are represented by horizontal lines, the left and right ends of which have 
the coordinates {\n.k +i /k Ti ,\n.k'^ i ) and (In — kn/k-i, ln/cf^) respectively. The phase space 
available for final-state emissions is given by area below the horizontal lines (including the 
fold that stick out of the main triangle). 

The ordering of the CCFM evolution in q + but not g_ means that this formalism is 
not left-right symmetric. In contrast, the LDC formulation is completely symmetric, which 
implies that the chain in Fig. 2 can be thought of as evolved from either the photon or the 
proton end. Thus, the LDC formalism automatically takes into account resolved virtual 
photon contributions. 

It is important to note that in the LDC model non-leading effects such as quark-initiated 
chains, the non-singular terms in the splitting functions and energy-momentum conservation 
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can be included in a straight-forward manner. The fact that fewer gluons are considered 
as initial-state radiation implies that typical z-values are smaller, thus resulting in smaller 
sub-leading corrections. It is also implies that in the LDC formalism the gluon distributions 
is quite insensitive to the second scale q, and to leading ml/a; it depends only on a single 
scale k?- Just as for CCFM, for finite x Sudakov form factors are introduced to regularize 
the 1/(1 — 2) poles in the splitting functions, thus introducing a q dependence. This de- 
pendence comes mainly from the last Sudakov form factor related to the convolution of the 
unintegrated parton densities with the off-shell matrix element, and the factorization 

$(x,4,g 2 ) = A s (g 2 )$(x,4) (9) 

is a good approximation for the LDC model. 



III. NUMERICAL RESULTS 



In this section we present the numerical results of our calculations and compare them 
with data from D0 and CDF. We will use here the expressions for the heavy quark hadropro- 
duction cross section and gluon-gluon fusion off-shell matrix element which were obtained 

inQ,Q- 

To illuminate the effect of the different contributions we will study the four different 
versions of LDC unintegrated gluon distributions presented inj^j, namely the standard, 
gluonic, gluonic-2 and leading ones. 

The standard version includes non-leading contributions from quarks and non-singular 
terms in the splitting functions. The gluonic and leading versions do not take into account 
quark links in the evolution. Furthermore the leading version does not includes non-singular 
terms in the splitting functions in contrast to the gluonic versions. For all versions, non- 
perturbative input parton densities of the formj^] 

xf i (x,l% Q )=A i x at (l-x) hi (10) 



where used, with all parameters Aj, Oj, b$ (i = g,d,d,u, . . .) and the perturbative cutoff k\ 



fitted 



TO 



to reproduce the measured data on i^x, Q ). The standard version has been fitted 
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in the region x < 0.3, Q > 1.5 GeV to experimental data taken by the HI [37], ZEUS 



;Q and E665Q 



NMC|3J| and E665 [40j collaborations. The gluonic and leading versions have been fitted in 
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FIG. 3: The prediction for the B meson px spectrum for \db\ < 1 at ^/s = 1800 GeV compared 
to the CDF data. The solid line is gluonic, dashed is glionic-2, short dashed is leading and dotted 



line is standard versions. Experimental data are from CDF 



28] 



the region x < 0.013 and Q 2 > 3.5 GeV 2 to data taken by the HI collaboration only. To 
study the sensitivity to the b parameter in equation (jlO)) an additional fit for gluonic case 
(called gluonic-2) was obtained (see [3] for detail information). We note that all versions 
give a satisfactory fit to the F 2 data. 

After we fixed the choice of the unintegrated gluon distribution, our theoretical results 
depend on the choice of heavy quark mass value, factorization scale /x 2 and selection of the 
heavy quark fragmentation function. In the present paper we convert heavy quarks into D 
and B mesons using the usual Peterson fragmentation function)^ with e = 0.06 for charm 
and e = 0.006 for bottom respectively^]. Also we used the following choice m c = 1.5 GeV, 



4.75 GeV and /x 2 



rrin 



m Q + Pt with Pt being the transverse momentum of the 



heavy quarks in the pp cm. frame. These choices are similar to the ones in e.g. 

The results of our calculations for bottom production are shown in figures 0HZI Figure 
El shows the prediction for the B + meson px spectrum for \ub\ < 1 at \fs = 1800 GeV 
compared to the CDF dataj^. One can see that results obtained with the leading version 
agree with the CDF data within experimental uncertainties. The gluonic version is very close 
to the leading one but goes a bit below the data. The results obtained using the gluonic-2 
and standard versions underestimate experimental data by a factor about 2 and are close 



to the NLO pQCD calculations 



291 ] . One can also see a difference in the shapes between 
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(a) 



(b) 
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gluonic — 
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D0 \-- 




FIG. 4: The cross section for muons with p?^ > 5 GeV (a) and p?^ > 8 GeV (b) from B meson 
decay as a function of rapidity compared to the D0 data. All curves are the same as figure 01 
Experimental data are from D0 j^l 

the standard version and the other ones, while the shapes of the gluonic and leading curves 
practically coincide. 



The recent D0 



31 



32j | experimental data refer also to muons which originate from 



the semileptonic decays of 5-mesons. To produce muons from B mesons in theoretical 
calculations, we simulate their semileptonic decay according to the standard electroweak 
theory. Our calculation of the cross section for muons from B meson decay as a function of 
rapidity da/d\y^\ is shown in figure |U for both p Tfl > 5 GeV (JI^) and pt^ > 8 GeV 
We find that only leading and gluonic versions agree with the D0 experimental data 
within errors. Again, the leading is somewhat above the gluonic although we now see that 
the difference due to the non-singular terms is mainly present in the central rapidity region. 
Also it is interesting to note that standard version have a more flat behavior compared 
to the other versions of the LDC unintegrated gluon distributions. The reason that the 
gluonic-2 results falls much faster at large rapidities is that it has a more rapidly falling 
input distribution at large x (b g = 7 in equation IjlOj) rather than b g = 4 as for the other 
versions') . Note that the NLO pQCD calculations underestimate data by a factor about 

Figure El shows the leading muon pt spectrum in the central rapidity region for bb produc- 
tion compared to the D0 data[31]. The leading muon in the event is defined as the muon 
with highest transverse momentum. Again we find that the leading and gluonic versions 
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FIG. 5: The leading muon pt spectrum for bb production compared to the D0 data. The cuts 
applied to both muons are 4 < pxa < 25 GeV, \n u \ < 0.8 and 6 < m uu < 35 GeV. All curves are 

n 

the same as in figure |31 Experimental data are from D0 |3J. 



FIG. 6: The double differential cross sections for muons from B meson decay with 2.4 < |g/J < 3.2 
as a function of Pt^- All curves are the same as in figure |21 Experimental data are from D0 j^ . 

agree well with the data while the other two are below. However, the shape of all curves 
practically coincide. 

Also the double differential cross sections da/dpT^dy^ in the forward rapidity region 
2.4 < < 3.2 are well described by the /c^-factorization approach with the leading and 
gluonic unintegrated gluon distributions (figure |UJ). Again there is only slight variations in 
the shape of the different cureves. The NLO pQCD calculations underestimate the D0 
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FIG. 7: Azimuthal correlations for muons pairs with 4 < pxa < 25 GeV, IryJ < 0.8 and 6 < m„„ < 

n 

35 GeV. All curves are the same as in in figure Experimental data are from D0 |3J|. 



data|32j by a factor of 4. 

n 

In a previous paper [27J it was found that investigations of bb correlations such as the az- 
imuthal opening angle between b and b quarks (or between their decay muons) is a powerful 
test for the different unintegrated gluon distributions. This is because these quantities are 
sensitive to the relative contributions of different production mechanisms to the total cross 
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In the naive collinear gluon-gluon fusion mechanism, the distri- 
bution over the azimuthal angle difference A<p b i must be a simple delta function S(A(f> b i — it). 
Taking into account non-vanishing initial gluon transverse momenta leads to the violation 
of this back-to-back quark production kinematics in the /^-factorization approach. 

The differential bb cross section da/dAeb^ for central muons is shown in figure For the 
overall cross section we se the same trend as in previous figures, where the leading is above 
the gluonic which in turn is above the standard. However, we here see from the shape that 
the leading clearly overestimates the decorrelation. The others curves agree better with the 
shape of the experimental result, although none of them are able to fully reproduce the peak 



at A0 



7T. 



Very recently the CDF collaboration have reported preliminary experimental datajii 
charm production at the Tevatron Run II. These data found to be about a factor of 1.5 
larger than NLO pQCD theoretical predictions Q]. data 3^ |. 

The results of our calculations of the transverse momentum spectra for centrally produced 
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FIG. 8: The prediction for the px spectrum of D* (a), D + (b) and D° (c) mesons with \vd\ < 1 at 
y/s = 1960 GeV compared to the CDF data. All curves are the same as in figure H3 Experimental 
data are from CDF[3J|. 

D*, D + and D° mesons at \fs = 1960 GeV are shown in figure |HJ Comparing with the B 
meson spectrum in figure H3 we see that the difference between the different versions of the 
LDC densities are the same both for the magnitude (i.e. leading > gluonic > gluonic-2 > 
standard) and for the shape (standard falls more steeply than the others). However, the 
data is now a bit closer to the calculations, especially for the leading version. 



IV. CONCLUSIONS 

We have presented results for charm and bottom production in pp collisions at high en- 
ergies within the framework of fey-factorization, using different unintegrated gluon densities 
obtained from the LDC model. It has previously been noted that the standard collinear 
factorization approach fails to describe the amount of heavy quarks produced at the Teva- 
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tron, and that /^-factorization may be the more correct way of describing the underlying 
physics. Our results agrees with this observation. However, within the LDC model it is also 
possible to study non-leading effects in the evolution, and it is clear that these introduce 
large uncertainties. We have found that the inclusion of quarks in the evolution has a big 
effect and that the results are sensitive to the treatment of non-singular terms in the gluon 
splitting function. Also we have found that the results are sensitive to the assumed shape 
of the non-perturbative input density at large x. 

In particular we find that we can only get an acceptable result for the case where only 
gluons are considered in the evolution and where only the singular terms in the gluon 
splitting function are included {leading). This is consistent with the findings in j^J where 
the CASCADE Monte Carlojs5] was shown to reproduce the amount of bottom production 
at the Tevatron. This program is also based on the CCFM equation and implements purely 
gluonic evolution without considering non-leading contributions in the splitting function. It 
is interesting to note that for the forward jet rates at HERA the situation is similar, as these 
can be reproduced by LDC and CASCADE only if non-singular terms are omitted jlfl|. 

Although the leading version of the LDC unintegrated gluon density gives a good de- 
scription of the overall rate of heavy quark production, less inclusive observables such as 
azimuthal correlations between bb (or pairs, are not as well described. 

Besides the uncertainties due to non-leading terms in the evolution, there is, of course, 
also some uncertainties due to the chosen values of the heavy quark masses. However, the 
effects of the heavy quark masses would not be large enough to change the conclusions 
presented here. In addition, there is some freedom in the choice of factorization scale. We 
have here chosen u = % which is the natural scale within the LDC model, but it can be 

n 

noted that in e.g. [25J, the scale was chosen to be /i = qr- Using such a scale in the LDC 
model would increase the bottom cross section, but, again, it would not be large enough to 
affect our conclusions. 

It is important to note that we have here only considered gluon-gluon fusion into heavy 
quarks. One may also expect non- negligible contribution from two other production mecha- 
nisms, namely gluon splitting, where a hard final-state gluon splits into a heavy quark pair 
in a subsequent final-state cascade, and heavy quark excitation processes, where the the 
hard process is gb — > gb, with the b produced by initial-state gluon splitting earlier on in 
the evolution. In principle, all these contributions can be taken into account in the LDCMC 
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Monte Carlo|36[ which implements the LDC model, however, this program is not yet able 
to fully handle hadron-hadron collisions 4s| . 

To conclude, it is most likely that fc^-factorization holds the key to understanding heavy 
quark production at the Tevatron. However, there are still large uncertainties, and much 
more work needs to be done before these processes are fully understood. 

V. ACKNOWLEDGMENTS 

The study was supported in part by RFBR grant N° 02-02-17513 and the Crawfoord 
Foundation (Sweden). A.L. was supported by INTAS grant YSF 2002 N° 399. 



[1] V. N. Gribov and L. N. Lipatov, Yad. Fiz. 15, 781 (1972). 
[2] L. N. Lipatov, Sov. J. Nucl. Phys. 20, 94 (1975). 
[3] G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977). 
[4] Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977). 

[5] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 44, 443 (1976). 

[6] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977). 

[7] I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 

[8] S. Catani, M. Ciafaloni, and F. Hautmann, Nucl. Phys. B366, 135 (1991). 

[9] J. C. Collins and R. K. Ellis, Nucl. Phys. B360, 3 (1991). 
[10] B. Andersson et al. (Small x) (2002), hep-ph/0204115. 
[11] M. Ciafaloni, Nucl. Phys. B296, 49 (1988). 

[12] S. Catani, F. Fiorani, and G. Marchesini, Phys. Lett. B234, 339 (1990). 

[13] S. Catani, F. Fiorani, and G. Marchesini, Nucl. Phys. B336, 18 (1990). 

[14] G. Marchesini, Nucl. Phys. B445, 49 (1995), hep-ph/9412327. 

[15] J. R. Forshaw and A. Sabio Vera, Phys. Lett. B440, 141 (1998), hep-ph/9806394. 

[16] B. R. Webber, Phys. Lett. B444, 81 (1998), hep-ph/9810286. 

[17] G. P. Salam, JHEP 03, 009 (1999), hep-ph/9902324. 

[18] B. Andersson, G. Gustafson, and J. Samuelsson, Nucl. Phys. B467, 443 (1996). 



15 



[19] B. Andersson, G. Gustafson, and H. Kharraziha, Phys. Rev. D57, 5543 (1998), hep- 
ph/9711403. 

[20] E. M. Levin, M. G. Ryskin, Y. M. Shabelski, and A. G. Shuvaev, Sov. J. Nucl. Phys. 53, 657 
(1991). 

[21] M. G. Ryskin, Y. M. Shabelski, and A. G. Shuvaev, Z. Phys. C69, 269 (1996), hep-ph/9506338. 

[22] M. G. Ryskin, A. G. Shuvaev, and Y. M. Shabelski, Phys. Atom. Nucl. 64, 1995 (2001), 
hep-ph/0007238. 

[23] P. Hagler, R. Kirschner, A. Schafer, L. Szymanowski, and O. Teryaev, Phys. Rev. D62, 071502 

(2000), hep-ph/0002077. 

[24] H. Jung, Phys. Rev. D65, 034015 (2002), hep-ph/01 10034. 

[25] S. P. Baranov and M. Smizanska, Phys. Rev. D62, 014012 (2000). 

[26] N. P. Zotov, A. V. Lipatov, and V. A. Saleev, Phys. Atom. Nucl. 66, 755 (2003). 

[27] S. P. Baranov, A. V. Lipatov, and N. P. Zotov (2003), hep-ph/0302171. 

[28] D. Acosta et al. (CDF), Phys. Rev. D66, 032002 (2002), hep-ex/0206019. 

[29] D. Acosta et al. (CDF), Phys. Rev. D65, 052005 (2002), hep-ph/0111359. 

[30] F. Abe et al. (CDF), Phys. Rev. D55, 2546 (1997). 

[31] B. Abbott et al. (DO), Phys. Rev. Lett. 84, 5478 (2000), hep-ex/9907029. 

[32] B. Abbott et al. (DO), Phys. Lett. B487, 264 (2000), hep-ex/9905024. 

[33] G. Gustafson, L. L6nnblad, and G. Miu, JHEP 09, 005 (2002), hep-ph/0206195. 

[34] C. Chen (2003), fERMILAB-THESIS-2003-14. 

[35] H. Jung and G. P. Salam, Eur. Phys. J. C19, 351 (2001), hep-ph/0012143. 

[36] H. Kharraziha and L. Lonnblad, JHEP 03, 006 (1998), hep-ph/9709424. 

[37] S. Aid et al. (HI), Nucl. Phys. B470, 3 (1996), hep-ex/9603004. 

[38] M. Derrick et al. (ZEUS), Z. Phys. C72, 399 (1996), hep-ex/9607002. 

[39] M. Arneodo et al. (New Muon), Phys. Lett. B364, 107 (1995), hep-ph/9509406. 

[40] M. R. Adams et al. (E665), Phys. Rev. D54, 3006 (1996). 

[41] C. Peterson, D. Schlatter, I. Schmitt, and P. M. Zerwas, Phys. Rev. D27, 105 (1983). 

[42] J. Chrin, Z. Phys. C36, 163 (1987). 

[43] M. Cacciari and P. Nason, Phys. Rev. Lett. 89, 122003 (2002), hep-ph/0204025. 

[44] M. Cacciari and P. Nason (2003), hep-ph/0306212. 

[45] G. Gustafson, L. Lonnblad, and G. Miu, Phys. Rev. D67, 034020 (2003), hep-ph/0209186. 



1(3 



